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► A model framework for entrained flow gasification is presented. 

► The model is validated using experimental data from various measurement setups. 

► The reaction regime of char conversion in industrial scale gasifiers is analysed. 

► At high temperature in the burner zone Regime III is approached. 

► In the later stages the reaction occurs under Regime II and Regime I conditions. 
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The paper presents a model framework that describes the reactions kinetics of entrained flow gasification. 
The model consists of different submodels for fuel devolatilisation, intrinsic char gasification, surface area 
evolution, char thermal annealing, pore diffusion, boundary layer diffusion, and particle size and density 
variation. These submodels are partly derived from measurements in previous publications and are partly 
based on theoretical derivations. The model framework is validated using fuel devolatilisation data and 
char gasification kinetics of a lignite that are measured in two entrained flow reactors at temperatures 
up to 1600 °C and pressures up to 1.0 MPa. A good correlation of model and experimental data is found 
for various operating conditions. 

The objective of the model development is the prediction of gasification behaviour in larger scale 
entrained flow gasifiers. To show its capabilities the model is applied to a one-dimensional plug flow reac¬ 
tor that represents a simple approach for a larger scale entrained flow reactor. Gas and temperature profiles 
within a 500 MW gasifier are predicted and cold gas efficiency and fuel conversion for different gasifier 
geometries and operating conditions are shown. 

Based on the model predictions the reaction regime of char conversion in entrained flow gasifiers is ana¬ 
lysed. At very high temperature in the burner zone, char gasification approaches Regime III conditions and 
a concentration gradient develops in the gas boundary layer around the char particles. The observed acti¬ 
vation energy is 62 kj/mol. In the medium stages of char conversion the reaction occurs under Regime II 
conditions and then approaches Regime I conditions at high conversion and at lower temperature. This 
is indicated by the high observed activation energy of 239 kj/mol in the final stages of conversion. 

As all three reaction regimes are relevant, fuel properties that have an impact on diffusion (e.g. particle 
size and pore structure) as well as on the intrinsic char reaction (e.g. intrinsic reactivity and thermal 
annealing) have to be considered to determine entrained flow gasification kinetics. 

© 2013 Elsevier Ltd. All rights reserved. 


1. Introduction 

Entrained flow gasification technology offers a wide range of 
applications for fossil as well as renewable energy sources. Solid 
fuels are converted at high temperature and high pressure and the 
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product gas can be used for power generation or as a raw material 
for chemical syntheses to produce liquid or gaseous fuels and chem¬ 
ical products. If C0 2 capture is required, a very efficient pre-combus¬ 
tion separation process is integrated. The possible decoupling of 
gasifier and final product output allows the installation of energy 
storage capacities that might be required due to the integration of 
fluctuating renewable energy sources into the electricity grid. 

Entrained flow gasifiers are the most important gasification 
technology. According to the 2010 Worldwide Gasification 
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Database [1] they account for most of the capacity installed 
worldwide. Solid fuels are converted at pressures of up to 
7.0 MPa and temperatures in the range of 1400-1600 °C by the 
addition of oxygen to the reaction zone [2], After devolatilisation, 
the volatiles are combusted in homogeneous gas reactions and 
heat is supplied for the subsequent endothermic heterogeneous 
char conversion. Char gasification processes with H 2 0 and C0 2 
are the slowest reactions and determine the rate of the overall 
conversion process [3-5]. 

Since char properties influence its conversion, it is important to 
understand typical char properties such as composition, pore 
structure, morphology, and surface area [6], The reactivity of char 
depends on total surface area, the number of reactive sites per unit 
surface area and the local gaseous reactant concentration [3], The 
conversion rate of char is influenced by the intrinsic reaction rate, 
but also by mass transport limitations within the char porous 
structure and the boundary layer. For a given coal char, the degree 
of limitation by chemical reaction or diffusion is mainly a function 
of reaction temperature. The dependence of the reaction rate on 
temperature is usually shown as a plot of reactivity (logarithmic) 
versus temperature (inverse), as shown in Fig. 1. This plot is gener¬ 
ally referred to as the Arrhenius plot and the influence of temper¬ 
ature is shown by the slope of the curve that gives the activation 
energy of reaction. 

The relation between chemical reaction and mass transport lim¬ 
itation is classified in three regimes. In Regime 1 at low temperature 
the chemical reaction rate at the char surface is very slow and deter¬ 
mines the overall reaction rate. The gas concentration is uniform 
within the char particle and the bulk gas phase. With increasing 
temperature the chemical reaction rate increases exponentially 
and becomes comparable to the pore diffusion rate within the par¬ 
ticle. With increasing reactant consumption at the inner particle 
surface a concentration gradient develops within the char particle. 
Under ideal Regime II conditions, the reactant concentration is 
equal to the bulk concentration at the outer particle surface and zero 
at the particle centre. When the temperature is increased further, 
the overall rate is controlled by the mass transfer between bulk 
phase and particle surface and a concentration gradient is formed 
in the boundary layer. As diffusion processes are only slightly influ¬ 
enced by temperature, the slope of the curve approaches zero. 

Measurements of reactions rates at low temperature under Re¬ 
gime I conditions are usually derived from thermogravimetric ana¬ 
lysers, and char conversion rates of different fuels are available 
(e.g. [7-10]). At higher temperature, the char gasification has to 



Fig. 1. Influence of the reaction temperature (inverse) on the char conversion rate 
(logarithmic) and on gas concentration gradients within a porous char particle. 


be measured in flow reactors under reaction conditions that are 
relevant to industrial scale entrained flow gasifiers. Only a few data 
sets are available that analyse char gasification at temperatures 
above 1000 °C. 

Ahn et al. [11] investigated the gasification rate of an Indone¬ 
sian sub-bituminous coal-char with C0 2 using a pressurised drop 
tube furnace (PDTF) at the Korea Electric Power Research Institute. 
Gasification rates were measured from 900 °C to 1400 °C and at 
1.0 MPa reactor pressure. From 900 °C to 1000 °C the activation en¬ 
ergy was determined to be 144 kj/mol. At higher temperature the 
activation energy decreased to 71.5 kj/mol which indicates the 
occurrence of mass transport limitations. The transition from Re¬ 
gime I to Regime II conditions occurred between 1000 °C and 
1100 °C. 

At the Central Research Institute of Electric Power Industry 
(CRIEPI), Japan, a pressurised drop tube furnace has been installed 
to analyse the gasification rate of coal char [12,13], The facility was 
operated at up to 1500 °C and at up to 2 MPa. For H 2 0 gasification 
of the NL coal char the transition temperature (Regime I-Regime II) 
is 1300-1400 °C, whereas it is 1200 °C for C0 2 gasification. For the 
S coal char no transition in C0 2 is observed up to 1200 °C. 

Matsumoto et al. [14] measured reactions rates of biomass char 
with C0 2 and H 2 0 in a pressurised drop tube furnace at the Naga¬ 
saki Research and Development Center, Japan. The PDTF was oper¬ 
ated at a total pressure of 0.4 MPa and at temperatures from 900 °C 
to 1300 °C. The measured activation energies were 136 kj/mol for 
the H 2 0 reaction and 94 kj/mol for the C0 2 gasification, both indi¬ 
cating the conversion under Regime II conditions. 

Coal gasification at a pressure of up to 2.0 MPa and tempera¬ 
tures of up to 1500 °C was measured in a pressurised entrained 
flow reactor at CSIRO [15,16], The activation energies measured 
at 0.5 MPa C0 2 partial pressure and a total pressure of 2 MPa were 
78-156 kj/mol. These were considerably lower than the activation 
energies of 242 kj/mol to 281 kj/mol that were measured for ex¬ 
tracted chars under Regime I conditions in a TGA and a fixed bed 
reactor. The transition temperature from Regime I to Regime II oc¬ 
curred around 1000-1100 °C dependent on char morphology. 

Other experimental gasification studies are carried out in atmo¬ 
spheric reactors [17-19], These measurements do not focus on the 
evaluation of the reaction regime and a transition temperature to 
Regime II is not given. 

The experimental data are limited to pressures of up to 2.0 MPa 
and temperatures of up to 1500 °C. A transition to Regime III is not 
reported under these conditions. Therefore, it is assumed that char 
conversion occurs under Regime II conditions in entrained flow 
gasifiers, as shown in Fig. 1 . Since the reaction temperature is usu¬ 
ally above 1500 °C on the industrial scale, especially in the burner 
zone of the gasifier, a transition to Regime III is possible and has to 
be considered when modelling entrained flow gasification. 

The objective of the paper is to measure gasification kinetics at 
high temperature and high pressure and extend the operation 
parameter range that is reported in the literature. The experimen¬ 
tal data are combined with devolatilisation and low-temperature 
char gasification measurements that are available from previous 
publications [4,5], Based on the experiments, a gasification model 
is derived and validated that describes entrained flow gasification 
in all three reaction regimes. Finally, the model is used to discuss 
the reaction behaviour within industrial scale gasifiers. 


2. Experimental methods and results 

2.1. Fuel preparation 

The proximate and ultimate analysis of the lignite (coal R) that 
is used in the analysis is shown in Table 1. The fuel is pulverised 















172 


1. Tremel, H. Spliethoff/Fuel 107 (2013) 170-182 


using a cross hammer mill and sieved between two meshes of 
width 80 pm and 160 pm using a tumbler screening machine. 
The moisture content is adjusted in a conditioning cabinet at a con¬ 
stant temperature of 40 °C and a defined humidity. After drying, 
the moisture content of coal R is 12 wt%. 

2.2. Entrained flow reactors 

Solid fuel gasification at high temperature is studied in two en¬ 
trained flow reactors. The Pressurised High Temperature Entrained 
Flow Reactor (PiTER) is described in detail in previous publications 
[4,5], The reactant gases (N 2 , H 2 , 0 2 , C0 2 , C 3 H 8 ) are supplied in gas 
cylinders and are fed to the top of the reactor by mass flow control¬ 
lers. Water steam is provided by a pressurised steam generator. 
Pulverised fuel is fed from a gravimetric dosing system with a 
highly controllable feed rate from 0.2 kg/h to 5.0 kg/h. The gas mix¬ 
ture is heated in the preheating tube and then enters the main 
reaction tube that has an inner diameter of 70 mm and a length 
of 2.3 m. The preheating tube and the reaction tube are surrounded 
by graphitic heating elements that allow a maximum operation 
temperature of 1800 °C. The reactor is situated within a pressure 
vessel that enables a maximum operation pressure of 5.0 MPa. 

The Baby High Temperature Entrained Flow Reactor (BabiTER) 
is designed for atmospheric pressure and on a smaller scale. A 
schematic diagram of this reactor is shown in Fig. 2. 

Different gases can be premixed and fed to the gas preheating 
unit. These gases are N 2 , 0 2 , H 2 , C0 2 , CO, CH 4 and air. The gas flow 
rate of each gas is adjusted by mass flow controllers. The gases are 
premixed in a gas collecting tube and then fed to the horizontal gas 
preheater. The ceramic preheating tube (inner diameter 25 mm, 
length 800 mm, 99.7% A1 2 0 3 ) is surrounded by a spiral-shaped 
electrical resistance heater that is connected to a transformer. 
The gas mixture can be heated to a maximum temperature of 
1300 °C and then enters the main reaction tube. 

A laminar flow is achieved by a flow straightener and the gas 
mixture streams downwards within the reaction tube (i.d. 
40 mm, length 1000 mm, 99.7% A1 2 0 3 ). The desired gas residence 
time within the reaction tube is achieved by adjusting the gas flow 
rate. The reaction tube is surrounded by two heating zones. Each 
heating zone consists of four MoSi 2 heating elements that are de¬ 
signed for a maximum temperature of 1800 °C. The temperature 
within each heating zone is controlled by a type B thermocouple 
that is situated at the outer wall of the reaction tube. 

Pulverised fuel is fed to the reaction tube from a dosing system 
at the top of the reactor. A constant fuel mass flow is achieved by a 
combination of screw feeder and vibrational chute and is control¬ 
lable in the range 50-500 g/h. The dosing system is situated on a 
balance allowing an online monitoring of the fuel feed rate. To pre¬ 
vent air from entering the reactor through the feed line, the feeding 
system is enclosed in a container. Inert gas fed to the container is 
used as a carrier gas for the fuel particles. The feeding system is 
connected to the reaction tube by a ceramic tube that ends after 
the flow straightener. 


Coal R 


Proximate analysis (wt.%, dry) 


Ultimate analysis (wt%, daf) 

Carbon 

Hydrogen 

Nitrogen 


67.8 

5.2 

0.8 

0.8 



The hot gas and particle mixture is directly cooled by spray 
water from a ring-shaped nozzle below the reaction tube. Char 
and gas samples are taken from the hot reaction tube through a 
water-tempered probe. This probe is inserted at the bottom of 
the quench zone and can be moved upwards into the reaction tube. 
Char and gas samples are taken at different stages of conversion by 
adjusting the height of the sampling probe. The samples are isoki- 
netically sucked in, whereby the gas flow rate is adjusted by a 
membrane pump. 

In both reactors an isothermal temperature profile is expected 
during the gasification experiments. In the PiTER experiments the 
oxygen concentration in the feed gas is only 2-7 mol%. The reac¬ 
tion of the coal particles with oxygen is exothermic which results 
in heat release. On the other hand, the inert gas in the reaction tube 
(nitrogen concentration is 93-98 mol%) can be seen as a thermal 
buffer which significantly reduces a possible temperature increase. 
Furthermore, the particle carrier gas and the fuel particles are in¬ 
jected into the hot reaction zone at a temperature of around 
25 °C (cooled injection probe). The temperature of the reaction 
tube at the point of fuel injection is monitored during the experi¬ 
ments. As a variation of temperature before and during coal injec¬ 
tion is not measured, a significant increase of temperature within 
the reaction tube is not expected and isothermal conditions are 
assumed. 

Due to the lower total pressure in the BabiTER experiments 
(atmospheric pressure), the oxygen concentration in the feed gas 
is higher (4-18 mol%) and consequently the thermal buffer capac¬ 
ity of the gas is lower. However, the reaction gas in the BabiTER can 
only be pre-heated to a maximum temperature of 1300 °C. In the 
experiments the pre-heating temperature is always chosen to be 
200 °C below the reaction tube temperature (an exception are 
experiments at 1600 °C where the maximum preheating tempera¬ 
ture of 1300 °C is used). As also the measurement of the reaction 
tube temperature at the point of fuel injection did not show an 
increase of temperature in gasification experiments, an isothermal 
temperature profile is also assumed in the BabiTER experiments. 

2.3. Determination of particle residence time 

The knowledge of the particle residence time is an important 
factor for the measurement of kinetic data, but cannot be directly 
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adjusted in an entrained flow reactor. The method for the determi¬ 
nation of the particle residence time in the experiments is briefly 
described in the following. 

The total gas flow rate in each experiment is used to adjust the 
average gas residence time in the reaction tube of both reactors. 
The average particle residence time can deviate from the gas resi¬ 
dence time. Possible deviations are caused by (i) the increase in gas 
amount by fuel gasification, (ii) incomplete entrainment of the par¬ 
ticles, and (iii) the flow velocity distribution within the reaction 
tube. 

The gas flow rate in the reaction tube continuously increases 
during the devolatilisation and gasification reactions. This leads 
to a decrease in the true gas residence time compared with the 
adjustments of the mass flow controllers. In a worst case calcula¬ 
tion, the devolatilisation occurs instantly at the point of fuel injec¬ 
tion and a higher gas flow rate exists along the whole reaction 
tube. In this case, the deviation between true average gas residence 
time and adjusted value is 0-7% for the PiTER experiments. In the 
atmospheric BabiTER experiments the influence of devolatilisation 
is higher. The true average gas residence time can be 2-22% lower 
than the set point residence time. Uncertainties in both reactors 
arise from the fact that it is not known a priori at which time 
and at which position in the reaction tube the increase in gas 
amount occurs. 

If the free fall velocity of the particles is significant compared to 
the gas velocity, the particles are moving faster than the surround¬ 
ing gas atmosphere. Both PiTER and BabiTER are designed as en¬ 
trained flow reactors where the gas flow rate is high. For each 
experimental run, Stoke’s law is applied and the terminal velocity 
of the fuel particles is calculated assuming a particle density of 1 g/ 
cm 3 , a diameter of 100 pm and a laminar flow regime. The entrain¬ 
ment is calculated from the deviation of absolute particle velocity 
and gas velocity. The particle entrainment in the PiTER and BabiT¬ 
ER experiments is between 82% and 97% with the average at 90%. 
This reduces the average particle residence by approx. 10% com¬ 
pared to the average gas residence time. 

The maximum reduction of particle residence time is calculated 
for each experimental run by determining the laminar free fall 
velocity of the particles. The true reduction is assumed to be half 
the maximum deviation, and the uncertainty to both sides is 50% 
of the maximum deviation. For example, if the maximum reduc¬ 
tion of the particle residence time due to incomplete entrainment 
is 0.2 s in an experiment, the true deviation is assumed to be 0.1 s 
and the uncertainty is ±0.1 s. 

The Reynolds number of the gas flow within the reaction tube is 
calculated for each experimental run. In PiTER experiments, values 
of the Reynolds number range from 800 to 4000 with the average 
being 1800. As the critical Reynolds number that indicates the 
transition between the laminar and the turbulent region is 2320, 
most experiments are nominally performed in the laminar region. 
However, in these experiments turbulence occurs during fuel 
injection and therefore in the upper part of the reaction tube. A 
laminar flow profile is only expected at the lower part of the reac¬ 
tion tube. As the gas velocity in the centre of a laminar flow profile 
is higher than the average flow velocity, particles that are in the 
centre of the reaction tube move faster than particles that are near 
the wall. In the turbulent region, the flow profile is described by a 
plug flow where the gas velocity is equal along the radial position. 
The particle velocity is then not dependent on the radial position. 

As the flow profile of all experiments can be described neither 
by complete laminar nor by complete turbulent flow, an exact ana¬ 
lytical calculation of the particle velocity as a function of radial po¬ 
sition is not possible. A solution to this problem is the positioning 
of the sampling probe. If the sampling probe were placed exactly in 
the centre of the reaction tube, the residence time of collected par¬ 
ticles would be shorter than the average gas residence time in a 


laminar flow but equal to the average gas residence time in a tur¬ 
bulent flow. The sampling probe is therefore placed at a defined 
position between the centre of the reaction tube and the wall. Then 
in a laminar flow the average gas residence time in the cross sec¬ 
tion of the sampling probe at this selected position approaches 
the average gas residence time in the cross section of the reaction 
tube. Therefore the residence time of particles collected with the 
probe also approaches the average gas residence time in the reac¬ 
tion tube. In a turbulent flow the position of the probe does not 
influence the average residence time of the particles collected. 
Therefore the residence time of particles collected at this probe po¬ 
sition equals the average gas residence time independent of the 
flow regime. 

During a set of measurements displacement is possible due to 
the vertical movement of the probe that might also influence the 
radial position. Hence a deviation in the experiments is expected. 
The distance between the sampling probe and the reaction tube 
wall is 2 mm, which is the maximum positive deviation. For the 
maximum negative deviation 5 mm are assumed. In a worst case 
consideration a pure laminar flow is assumed which results in an 
uncertainty of the average gas residence time within the probe 
cross section of -33% and +8%. 

In BabiTER experiments the Reynolds number is between 90 
and 370. A laminar flow profile is expected in all experiments. 
Due to the flow straightener and the particle injection in flow 
direction, only small turbulences are likely in the upper part of 
the reaction tube. 

Because of the small fuel mass flow, the sampling probe is 
placed in the centre of the reaction tube to collect a sufficient char 
mass in an appropriate time interval. As the true gas residence 
time in the centre of a laminar flow profile is lower than the aver¬ 
age gas residence time in the whole cross section, the residence 
time of the collected particles is reduced. Due to the laminar flow 
in each experiment which yields identical flow profiles, the devia¬ 
tion between the particle residence time and the average gas resi¬ 
dence time is constant. The particle residence time is 41% lower 
than the average gas residence time. Due to an uncertainty in isoki¬ 
netic conditions that is caused by the imprecise measurement of 
the gas volume flow through the sampling probe, the uncertainty 
of the particle residence time is 8% for each experiment. 

Due to the combined effects of devolatilisation, incomplete 
entrainment, and the flow profile in the reaction tube, the true par¬ 
ticle residence time of the char collected differs from the average 
gas residence time that is adjusted in each run. The true particle 
residence time and the uncertainty are calculated for each experi¬ 
mental run. The uncertainties of the three effects are combined 
according to the classical Gaussian error propagation method [20], 

2.4. Measurement of char conversion 

The char samples collected during each experimental run of the 
entrained flow reactors are used to calculate elemental and overall 
conversion by an ash tracer method. It is assumed that the amount 
of solid mineral matter does not change during conversion. The po¬ 
tential devolatilisation of ash components at high temperature is 
neglected. The ash and elemental (CHNS) content of each char 
samples is measured and the conversion is calculated by a compar¬ 
ison to the initial fuel composition. The uncertainty in the analysis 
is determined by repeated measurements using coal R, a reference 
char (produced from coal R under the conditions of the proximate 
analysis) and a PiTER char sample (pyrolysis char of coal R, 1200 °C, 
0.5 MPa, 1.3 s residence time). In the coal analysis the uncertainty 
in ash content is ±0.14 wt% and in carbon content ±0.30 wt%. The 
char analysis reveals an uncertainty in ash content of ±0.15 wt% 
and in carbon content ±0.41 wt%. The uncertainties from the labo¬ 
ratory analysis are used to calculate uncertainties in the conversion 
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data. Following the Gaussian error propagation procedure [20] the 
uncertainty u of an output quantity y (e.g. conversion) is related to 
the uncertainties of input quantities Xk (e.g. char analysis): 



The confidence interval is defined by the uncertainty as follows: 


experiment with the char properties after the first experiments. 
As the char mass collected during PiTER experiments is very small 
(small sampling probe, isokinetical sampling conditions) compared 
with the required fuel mass for an experiment, experimental runs 
using char as the fuel are not possible with the current setup. 

2.5. Experimental results 


y - eu(y) < y < y + eu(y) 


(2) 


where e u(y) is the extended uncertainty. When a value of 2 is used 
for e, the real value of y lies in the confidence interval with a confi¬ 
dence level of 95%. The e = 2 criterion is often used in industrial 
metrology [20], Based on the measurement uncertainty of the lab¬ 
oratory analysis, the absolute extended uncertainty of carbon con¬ 
version and overall conversion is 2 wt% (95% confidence interval). 

The char sample mass after the BabiTER experiments is not al¬ 
ways large enough to measure ash content. The oxygen concentra¬ 
tion in a char sample after devolatilisation is very low. Therefore 
the ash content can be calculated directly from the CHNS analysis 
assuming zero oxygen concentration. This results in a higher 
uncertainty of the ash content and the average absolute extended 
uncertainty of carbon and overall conversion is 3 wt% (95% confi¬ 
dence interval). 

As pulverised coal is used in all gasification experiments, the 
measured conversion is a combination of devolatilisation and char 
burn-off. The conversion of char cannot be measured directly in 
these experiments. Therefore pyrolysis (N 2 atmosphere) and gasi¬ 
fication (O2/N2 atmosphere) experiments are carried out for each 
set of operation conditions. The char conversion x char in the gasifi¬ 
cation experiment is then calculated by the subtraction of the car¬ 
bon conversion in the corresponding pyrolysis experiment: 


Xchar 


Xcgas Xc.pyr 


(3) 


where xc, gas is the carbon conversion in the gasification experiment 
and x CiP yr is the carbon conversion after pyrolysis in a nitrogen 
atmosphere. 

Alternatively to the experiments shown here, char conversion 
could be measured directly if the reactor is operated with prepared 
char. This would require a first experiment for the production of 
char under defined conditions (pyrolysis experiment). During this 
experiment char is collected and is used as fuel in a second exper¬ 
iment. A direct determination of char conversion in the second 
run is possible by comparing the char analysis after the second 


Gasification experiments are carried out in both entrained flow 
reactors. In all experiments a constant stoichiometry is applied. For 
each carbon atom in the fuel mass flow one oxygen atom in the 
reaction gas is provided. This results in a constant oxygen to carbon 
ratio (O/C) of 1.0 ±0.1 The uncertainty results from deviations of 
the gas and the fuel mass flow rate. Nitrogen is used as the balance 
gas and as fuel carrier gas. The PiTER is operated with a mass flow 
rate of 1.25 kg/h and during the BabiTER operation the coal flow 
rate is 0.2 kg/h. The temperature in both reactors is varied between 
1200 °C and 1600 °C and the PiTER is operated at 0.5 MPa and 
1.0 MPa. 

The entrained flow gasification data of coal R are shown in 
Fig. 3. Experiments at a low particle residence time and under 
atmospheric pressure are carried out in the BabiTER. In the PiTER 
the residence time is increased to 2.4 s. 

Both overall and char conversion increase rapidly with particle 
residence time. The temperature has a significant influence on 
the conversion. However, the difference in conversion between 
1200 °C and 1400 °C is small and decreases further between 
1400 °C and 1600 °C. In the early stages of conversion there is a 
steep increase in char conversion indicating a fast reaction rate. 
After a residence time of 0.5 s about 50% of the char is already gas¬ 
ified. The reaction rate slows down at higher residence times and 
conversions. 

There is a continuous trend in the data from 0.1 MPa to 1.0 MPa 
when particle residence time is increased. This is due to the fact 
that the data for a constant temperature lie within a small band. 
This trend is independent of the operation pressure indicating no 
significant influence of the total pressure on reaction rate. 


3. Modelling approach and validation 

The reaction behaviour of coal R has been analysed in more de¬ 
tail in previous publications. Fuel devolatilisation and the loss of 
char reactivity at high temperature and longer time at temperature 
have been analysed under entrained flow conditions [4], 


(a) 



0 0.5 1 1.5 2 2.5 

Residence time [s] 


0 0.5 1 1.5 2 2.5 

Residence time [s] 


O 1200°C; 0.1 MPa 

♦ 1200°C; 0.5 MPa 
X 1200°C; 1.0 MPa 
A 1400°C; 0.1 MPa 
A 1400°C; 0.5 MPa 
O 1600°C; 0.1 MPa 

• 1600°C; 0.5 MPa 


Fig. 3. Conversion of coal R during entrained flow gasification in the PiTER (0.5 MPa and 1.0 MPa) and in the BabiTER (0.1 MPa) at different temperature and particle residence 
time: (a) overall conversion; (b) char conversion (error bars are shown exemplary for the 1200 °C data and denote the uncertainty in the experiment). 
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A devolatilisation model and a char deactivation model are derived 
based on experimental data. The influence of operation conditions 
on the heterogeneous char gasification under Regime 1 conditions 
is analysed in Ref. [5], Furthermore, the change of char surface area 
is analysed under entrained flow gasification conditions [5], Mod¬ 
els for the intrinsic char reaction rate and the surface area develop¬ 
ment are available. First, the submodels derived in previous 
experiments are briefly described and then the modelling approach 
is further extended with submodels for pore diffusion, boundary 
layer diffusion and the variation in particle size and density during 
conversion. All model parameters are given in Table 2. 

3.1. Description of the submodels 


3.1.1. Devolatilisation submodel 

The influence of pressure and temperature on volatile yield dur¬ 
ing devolatilisation is described by [4] 


Yv(p) = Y v , Pstt 


In (P/Pset) 
<5 


(4) 


Yv(T) = Yy,r sa + (Y v , Tma - Y VJset ) ■ (1 - exp[-0(T - T«t)]) (5) 

where 5 is a measure for the influence of pressure and p set is a ref¬ 
erence pressure. The parameter i? is a measure for the influence of 
temperature, T set is a reference temperature, and T max is the maxi¬ 
mum temperature where the maximum volatile yield Yyjmax is 
achieved in the experiment. 

The rate of devolatilisation is described by a single first order 
rate equation 

dY jf 1 = Ay ■ exp (Y vfmal - Y v (t)) (6) 


different submodels. 


Pyrolysis submodel 


0.5 MPa 
1000°C 
0.564 


Intrinsic reactivity submodel 


(475 + 37 p[MPaJ) m 2 /g 


Char deactivation submodel 
Amax 23.4 

F„ 21.0 X to 3 s-’ 

E af 117.2 kj/mol 

Particle size & density submodel 
Swelling index 1 

Fuel particle diameter 50 pm 

Initial char density 1.04 g/cm 3 

Pore diffusion 
Molecular diffusivity 

Tortuosity x 
Pore diameter 




Measurement 

Measurement 

Measurement 


Measurement 


where Ay is the pre-exponential factor, E AV is the activation energy, 
and Yv final is the final volatile yield at operating pressure and tem¬ 
perature given by (Eqs. (4) and (5)). 

3.1.2. Intrinsic reactivity submodel 

The intrinsic reactivity at the char surface is modelled by an nth 
order rate equation 

r intr = k 0 • exp -p? (7) 

where the rate is described by the pre-exponential factor k 0} the 
activation energy E A , and the exponent n [5], The equation is used 
to describe the char reactions with C0 2 and H 2 0 where Pi is the par¬ 
tial pressure of the reactant gas. 

3.1.3. Char surface area submodel 

The initial char surface area after devolatilisation is described 
by an empirical linear approximation of experimental data at dif¬ 
ferent pressure [5] 

S jn , [m 2 /g] =475 + 37 p [MPa] (8) 

The change of surface area with conversion is based on the Ran¬ 
dom Pore Model (RPM) [21,22], However, the experimental surface 
area data from the PiTER do not correlate with the mass specific 
form of the RPM [5], However, when the mathematical description 
of the volumetric surface area [m 2 /m 3 ] is used to predict mass spe¬ 
cific surface area [m 2 /g] a good correlation with experimental data 
is found. 

The surface area during conversion is modelled by: 

s Ff] = s,ni ■ ° ■ x) ■ V 1 - ? ' in ( i -*) « 

where the initial surface area S, n i is given by (Eq. (8)) and 'P is the 
only adjustable parameter. 

3.1.4. Char deactivation submodel 

A significant decrease in the intrinsic reactivity is detected dur¬ 
ing heat treatment and char conversion [4], In the experiments it 
was observed that thermal deactivation does not significantly 
influence the activation energy of the intrinsic char conversion, 
but has a major effect on the pre-exponential factor [5], Thermal 
deactivation is dependent on time t and temperature T and is con¬ 
sidered in the model by an extension of the pre-exponential factor 
in (Eq. (7)): 


fco = [1 +/ • (Amax - 3)] • kgp 

(10) 

/ = exp[-F„.exp(^)t] 

(11) 


The parameter A max is a fuel-specific constant. The parameter 
k 0 .D is the pre-exponential factor of a strongly deactivated char 
sample in the intrinsic model and is derived from experimental 
data. 

The pre-exponential factor F 0 and the activation energy E AF de¬ 
scribe the loss of active sites by a first order Arrhenius type equa¬ 
tion [4]. 

3.1.5. Pore diffusion submodel 

During char gasification at high temperature, chemical surface 
reactions occur simultaneously to mass transfer within the porous 
char structure. Due to mass transfer limitations a gas concentration 
gradient within the char particle might be established and the sur¬ 
face in the particle centre is exposed to a lower reactant concentra¬ 
tion than the surface near the particle exterior. As the reaction 
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within a porous structure is comparable to heterogeneous cataly¬ 
sis, most of the theoretical development is similar. 

The degree of mass transport limitation is described by the 
effectiveness factor r\ that is defined as the ratio of observed reac¬ 
tion rate and proposed reaction rate if the reactant concentration 
at the particle outer surface is constant throughout the particle 
[23], 

The effectiveness factor is a function of the Thiele modulus 0 
and can be described for a spherical particle shape [24] 

^ = 0 [tanh(30) ~ 3</>] (12) 

Bischoff [23] developed a general expression for the Thiele 
modulus in a semi-infinite particle, referred to as the general 
modulus 


surface concentration is linked to bulk gas concentration by mass 
transfer in the boundary layer of a char particle. This second 
diffusion limitation is not considered by the effectiveness factor 
approach. Therefore a further submodel is required to account 
for mass transport limitation within the boundary layer of a char 
particle. 

This model is based on a simple form of Fields Law [3]: 

J i = hM i (c bm -c s ,i) (15) 

where J is the mass flux of reactant gas i at the particle outer sur¬ 
face, h is the mass transfer coefficient, M,- is the molar mass of the 
reactant gas, and c buUCii and c Sfi are molar concentrations in the bulk 
phase and at the outer particle surface, respectively. In order to re¬ 
late the gas diffusion rate to the reaction rate of char, a gravimetric 
stoichiometric coefficient N is introduced [3]: 




(13) 


(16) 


where L is a characteristic length, r is an arbitrary rate expression, c s 
is the gas concentration at the outer surface, D eff is the effective dif¬ 
fusion coefficient and a is a substitute for the integration. This def¬ 
inition will result in all solutions for the effectiveness factor 
coinciding for large <P and the differences are relatively small over 
the whole range of <P [23]. Earlier, Aris [24] used particle volume 
V P and external surface area S P to express the characteristic length 
L = Vp/Sp. 

The Thiele modulus for the char conversion is derived by com¬ 
bining (Eqs. (7) and (13)). The gas concentration Cs is substituted by 
gas partial pressure and the following equation is derived: 

where d is the char particle diameter, p is the apparent char density, 
S is the mass specific char surface, D eS f is the effective diffusivity and 
Ps,i is the reactant gas partial pressure at the char outer surface. The 
effective diffusivity is derived from molecular diffusivity, Knudsen 
diffusivity, particle porosity e, and pore tortuosity factor x. The for¬ 
mulas for the calculation of effective diffusivity are taken from Ref. 
[25], 

The initial particle porosity after devolatilisation is calculated 
from the measured particle density of a pyrolysis char and the den¬ 
sity of graphite (2.27 g/cm 3 [26]). A typical value of 3 is taken for 
the tortuosity factor [25], 

The effectiveness factor is then calculated by using (Eq. (12)). 
This equation is initially derived for a first order reaction. The 
effectiveness factor for the nth order reaction can be approximated 
by the first order curve [3], Hong [27] calculated the maximum er¬ 
ror that occurs from the approximation by the comparison to an 
exact numerical solution. The maximum error is 16% for <P = 0.7 
and n = 0, but is negligible for 0 < 0.2 and 0 > 5. For n > 0.25 the 
maximum error is always less than 10% which is valid for the val¬ 
ues of n used here. 

The researchers from CSIRO [16,28,29] derived similar formula¬ 
tions to account for pore diffusion and further included modifica¬ 
tions to account for the effect of convective flow in the pore 
structure. Convective flow only increases the general Thiele modu¬ 
lus if the diffusion rate is limited by molecular diffusivity which 
occurs at larger pore diameters. Since char has a microporous 
structure, Knudsen diffusion is more important and the effect of 
convective flow is not considered here. 

3.1.6. Boundary layer diffusion submodel 

The effectiveness factor approach and the calculation of reac¬ 
tion rate in the pore diffusion regime are based on the knowledge 
of the gaseous concentration at the particle outer surface. The 


where M c and M,- are the molar weights of carbon and the reactant 
gas, and v is the stoichiometric coefficient of the gasification reac¬ 
tion. N for the char-C0 2 reaction is 3/11, and for the char-H 2 0 reac¬ 
tion 2/3. The reaction rate of carbon at the outer particle surface is 
then 

r c,i,e« (17) 

The mass transfer coefficient h is obtained from a Sherwood 
Number correlation for a sphere [30]: 

Sh = ^= 2 + 0.6Re° 5 Sc ° 33 (18) 

l>m, i 

In entrained flow gasification, gas and particles move at almost 
the same velocity and hence a stagnant atmosphere can be as¬ 
sumed around a char particle, i.e. Re a; 0 and Sh ~ 2. In this case, 
(Eqs. (17) and (18)) are used and gas partial pressure is introduced 
to calculate the observed reaction rate 


d 2 pRT 


- {Pbulk.i ~ Ps,i) 


(19) 


where d is the particle diameter, p the apparent density of the par¬ 
ticle, and Dm,! is the molecular diffusion coefficient for the gas i. The 
diffusion coefficient can be calculated for binary gas mixtures [25], 
Under ideal Regime III conditions, the intrinsic reaction rate is 
so fast that the partial pressure at the outer surface approaches 
zero. The observed reaction rate is then only dependent on bound¬ 
ary layer diffusion. 


3.1.7. Particle size and density submodel 

The change of particle size during devolatilisation is modelled 
by a swelling index that is the ratio of char particle diameter after 
devolatilisation and initial fuel particle diameter. For the coal R no 
significant change in particle size during devolatilisation is mea¬ 
sured and a swelling ratio of unity is used. 

Particle diameter d and density p can change during char con¬ 
version depending on the reaction regime. As particles only react 
at the outer surface under Regime III conditions, the particle size 
decreases but the particle density remains constant. In contrast, 
in Regime 1 the reaction occurs at the internal particle surface 
and the particle density decreases whereas the diameter almost 
remains constant. In coal combustion modelling power law rela¬ 
tions are often used to describe variations of the particle diameter 
and density with conversion [31]: 
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The subscript 0 denotes initial char conditions after devolatili¬ 
sation. For spherical particles, a.+ 3/1= 1, and the parameters are 
assumed to be constant during conversion [31], In Regime III, 
a = 0 and p = 1/3, and particles burn at constant density. In Regime 
I, a = 1 and p = 0, and the diameter remains constant. In entrained 
flow gasifiers, char conversion can occur in all three reaction re¬ 
gimes, thus 0 < a < 1 and 0 < p < 1/3. 

When neglecting the mineral matter content, the mass ratio can 
be described by char conversion: 


Then, the particle diameter and density are correlated by the ra¬ 
tio of a//?: 



Essenhigh [32] analysed an experimental combustion data base 
containing 37 values obtained from 15 literature sources. In spite 
of the different experimental methods, char types, fuel origins, par¬ 
ticle sizes and reaction temperatures, the ratio a IP was always in 
the range 0-7, and with most in the range 1-3. 

Liu and Niksa [33] modelled char gasification in entrained flow 
and fluidised bed reactors and used a constant value of 0.95 for the 
parameter a assuming that the reaction occurs almost completely 
under Regime I conditions. Monson et al. [34] measured char com¬ 
bustion at elevated pressure and particle temperatures from 
1123 °C to 1723 °C. The change of particle density with burn-off 
was best described by values for a between 0.5 and 0.65 indicating 
Regime II conditions. 

In Regime I the ratio a IP approaches infinity and under Regime 
III conditions it is zero. The reaction regime can be expressed by 
the quantity of internal and external char surfaces that are in¬ 
volved in the reaction. The internal surface area S int is larger by 
far than the external surface area S ex . In Regime II, the percentage 
of internal surface that is involved in the reaction is described by 
the effectiveness factor rj. The ratio of involved internal surface 
area and external surface area is used to describe oc/y5: 


a. 

J 


nSim 
3S m 


(24) 


In Regime I, rj = 1 and ot/yS approaches a large number. In Regime 
III, r\ S int ~ S c . xt and a.IP becomes 1/3. Thus the severity of pore dif¬ 
fusion in Regime II is described by the magnitude of r| and the ratio 
alp can be described in all three regimes. For a spherical particle, p 
is a function of a. and the reaction conditions: 


tySdp 0 (l-x)“ 


(25) 


where S is the specific surface area and d is the particle diameter. 

When using a constant value of a (for example 0.95 that is pro¬ 
posed by Liu and Niksa [33]), the parameter p can be calculated 
from the effectiveness factor and particle properties. 


3.2. Model validation using entrained flow data 


Due to the combination of several submodels and the imple¬ 
mentation of process parameters that are not constant during con¬ 
version (e.g. gas concentrations) numerical methods are used to 
combine the submodels and to calculate fuel conversion. 

The PiTER is assumed to be a one-dimensional plug flow reactor 
and is divided in cells in gas flow direction. For each cell, the fuel 
conversion rate is calculated based on process parameters. At the 


cell outlet the main gas concentrations (H 2 , CO, H 2 0, C0 2 , N 2 , and 
0 2 ) are calculated assuming chemical equilibrium. 

The fuel conversion is subsequently divided into three pro¬ 
cesses that occur consecutively in the reaction tube. In the first cell, 
gas and fuel particles are mixed. Gas concentrations are calculated 
from the inlet gas concentrations and from the fuel moisture con¬ 
tent. The devolatilisation process is modelled in the subsequent 
cells. In each cell the devolatilisation kinetics of the pyrolysis mod¬ 
el are applied and the total volatile yield after each cell is calcu¬ 
lated. The elemental composition of volatiles is based on the 
ultimate analysis (carbon, hydrogen, nitrogen and oxygen content) 
of the fuel. Sulphur, chlorine and other minor components are not 
considered. 

The focus of the model development is on the heterogeneous 
char reactions. After devolatilisation, char conversion occurs in 
the remaining length of the reaction tube. In each cell char conver¬ 
sion rate can be calculated using the instantaneous process param¬ 
eters (e.g. temperature, gas concentrations, pressure), char 
properties (e.g. surface reactivity, density, particle diameter) and 
chemical/physical properties (e.g. diffusivity). 

The only parameter in the model framework that cannot be 
derived from the measurement data is the effective pore diame¬ 
ter that affects the extent of pore diffusion limitation. The aver¬ 
age pore diameter in the simulation is used as an adjustable 
parameter to minimise the deviation between the experimental 
data and the model prediction. This procedure leads to an aver¬ 
age pore diameter of 47.3 A. The gasification reactions occur 
mainly within the micropores [35] that are defined by radii less 
than 10 A. However, during char conversion, the reaction will 
also occur in the mesopores (up to 200 A) and the average pore 
size will be shifted to larger diameters. The contribution of lar¬ 
ger pores to the total pore volume depends on carbon type 
[36] and conversion [35], A certain part of the pore volume 
can be in the range of 10-50 A [35] or even in pores with diam¬ 
eters up to 350 A [36], Therefore the value for the average pore 
diameter is in the expected range. 

The comparison of simulation and experimental data after 
pore size adjustment is shown in Fig. 4. The deviation between 
simulation and experimental data is small and the fuel conver¬ 
sion can be predicted by the numerical calculation based on 
the different submodels. As the simulation predicts the complete 
gasification process and considers changes in the gas phase as 
well as particle properties, further simulation results can be 
compared with experimental data. Fig. 5a compares the model 
predictions of the C0 2 yield (mol C0 2 per mol initial carbon sup¬ 
plied by the coal) in the PiTER reaction tube to the C0 2 yields in 
the experiments. The experimental C0 2 yield is calculated from 
C0 2 concentrations measured by mass spectroscopy during 
each run and from the total gas flow rate within the reaction 
tube. 

There is a good correlation of the C0 2 yields between simulation 
and experiment. Pyrolysis and volatiles combustion lead to a max¬ 
imum in C0 2 yield at short residence time. In the subsequent het¬ 
erogeneous reactions, C0 2 is consumed and the yield continuously 
decreases. As the fuel conversion increases at higher temperature, 
the C02 yields are decreasing from 1200 °C to 1600 °C. 

The evolution of specific surface area during char gasification 
is shown in Fig. 5b. The measurement uncertainty of the surface 
area (C0 2 -DFT) is low and there is a very good predictability of 
the simulation. The small increase in surface area that is ob¬ 
served in BabiTER experiments and also expected in the PiTER 
at shorter residence time is predicted by the simulation. At long¬ 
er residence times (higher conversion) the surface area continu¬ 
ously decreases. The deviation of the 1600°C data is caused by 
the overestimation of char conversion in the simulation that is 
shown in Fig. 4a. 




Fig. 4. Comparison of the model prediction and experimental data for the gasification of coal R in the PiTER: (a) char conversion; (b) overall conversion (the BabiTER data at 
atmospheric pressure are shown for comparison). 



Fig. 5. Modelling of coal R gasification in the PiTER at 0.5 MPa compared to experimental data (BabiTER data at atmospheric pressure are shown for comparison): (a) C0 2 
yield in mol C0 2 per mol of carbon; (b) specific surface area. 


4. Design and optimisation of larger scale entrained flow 
gasifiers 

4.J. Description of the simulation approach 

A major difference between the PiTER facility and a large scale 
gasifier is the heating of the reaction zone. During the PiTER oper¬ 
ation the reaction tube is kept at a constant temperature to observe 
reaction rate at isothermal conditions. The implementation of an 
energy balance for each cell is required to extend the model to lar¬ 
ger scale gasifiers. For this purpose, the reaction enthalpy in each 
cell is calculated and the local gasifier temperature is derived using 
heat capacities of gas and solids. Values of the heat capacities of 
gases at reaction temperature are taken from the NIST database 
[37], The heat capacity of char is approached by the heat capacity 
of graphite [38] and the heat capacity of ash is estimated to be 
1000J/(kgK). In the devolatilisation stage the heat capacities are 
taken at a temperature of 1500 K. Standard enthalpies of formation 
of H 2 0, C0 2 and CO [39] are used to calculate the reaction enthalpy. 


The enthalpy of vaporisation of water is 2256 kj/kg [40] and the 
ash melting enthalpy is assumed to be 1000 kj/kg. The temperature 
is first increased by the exothermic volatiles combustion and then 
continuously decreases due to the endothermic char-H 2 0/C0 2 
reactions. Heat losses through the gasifier wall are not considered 
and an adiabatic operation is assumed. 

Due to the usage of a plug flow model and further simplifica¬ 
tions there are some limitations to the simulation. The one-dimen¬ 
sional plug flow model does not consider back mixing effects and 
turbulence. As particle buoyancy and inertia are not considered, 
the particles are assumed to be completely entrained in the gas 
flow. The assumption of chemical equilibrium in the gas phase at 
high temperature is a good approximation. However, only the main 
gas components are considered. The very high temperatures after 
the volatiles combustion enable further reactions like molecule 
separations, atomisation and ionisation that would lower the en¬ 
thalpy of reaction and also the maximum temperature, but are 
not considered in the simulation. During char conversion, the frag¬ 
mentation of particles is expected. For lack of experimental data, 
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fragmentation is not considered in the simulation. Ash reactions 
are only considered by an ash melting enthalpy and a heat capac¬ 
ity. Chemical reactions of mineral compounds or the limitations of 
the gasifier operation by slag properties are not considered. 

The limitations result in an overestimation of reaction temper¬ 
ature in the flame zone. Due to the absence of turbulence and 
the non-consideration of high temperature gas reactions, the 
temperature in the combustion zone and in the first stages of char 
conversion is overestimated. As the reaction rates are temperature- 
dependent, an over-prediction of fuel conversion is expected. 
However, the char conversion approaches Regime III conditions 
at the highest temperatures where the temperature has only a 
minor influence on reaction rate (see Fig. 1 ). Therefore, the overes¬ 
timation of temperature does not significantly increase the char 
conversion. This has been proved by several test calculations. In 
the medium and later stages of conversion, a good predictability 
of temperature is expected. 
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Fig. 6. Influence of gasifier geometry (diameter) on the operation performance 
(cold gas efficiency and fuel conversion) and outlet temperature of a 500 MW 
entrained flow gasifier for the conversion of coal R (fixed gasifier length to diameter 
ratio of 1.8). 


4.2. Basic design of an entrained flow gasifier 

The gasifier diameter and its length are the size parameters in 
the simulation. Both can be adapted to the thermal fuel input, to 
the operation pressure and to the fuel properties. The design aims 
at complete fuel conversion and minimisation of the reactor size. 
The gasifier has an optimum cold gas efficiency and the ash/slag 
is free from organic matter only if the fuel conversion is high. 
The input parameters for the simulation are taken from Table 2 
and gasifier operation parameters are shown in Table 3. The aver¬ 
age fuel particle size is set to 100 pm. 

Schingnitz and Mehlhose [41 ] give basic geometries of the GSP 
entrained flow gasifier. Diameter and length increase with increas¬ 
ing thermal fuel input, but the ratio of length to diameter is almost 
constant at 1.8. 

For the selection of an optimum gasifier size for a thermal fuel 
input of 500 MW the gasifier diameter is increased and the length 
is adapted. The stoichiometry is defined by an O/C ratio of 0.95 
which corresponds to an oxygen mass flow rate of 65.8 t/h (fuel 
feed rate 90.3 t/h). The cold gas efficiency, fuel conversion, and gas¬ 
ifier outlet temperature for different geometries are shown in 
Fig. 6. 

A reduction of the gasifier size reduces the reaction time and 
therefore the fuel conversion. The cold gas efficiency is hardly 
influenced down to a gasifier diameter of 2.25 m (length 4.05 m, 
residence time 1.45 s). A further decrease in the diameter results 
in considerable losses due to incomplete conversion. Since a lower 
amount of thermal energy is required for the endothermic char 
reaction if the conversion is incomplete and since heat losses are 
not considered, the gasifier outlet temperature increases when 
the diameter is decreased. 

The gasifier diameter of 2.25 m results in a high fuel conversion 
of 99.6%. The gas concentrations and temperature profile within 
the first meter of the gasifier are shown in Fig. 7. The change of 
gas concentrations and temperature is smaller in the remaining 
length of the gasifier. 


Table 3 

Operation parameters for the simulation of an industrial scale entrained flow gasifier. 



Fig. 7. Gas concentration, conversion and gas temperature within an industrial 
scale entrained flow gasifier (gasifier length 4.05 m, diameter 2.25 m). 


The temperature after the volatiles combustion reaches a max¬ 
imum of 2430 °C which might be slightly overestimated, but then 
decreases rapidly. Despite the mass transport limitations, the reac¬ 
tion rate is initially very high and most of the fuel and char conver¬ 
sion occurs in the first metre of the gasifier length. Thermal 
annealing of char is very fast at the highest temperatures and the 
char is completely deactivated after about 0.5 s of the total resi¬ 
dence time of 1.45 s. The thermal deactivation, the loss of surface 
area, and the continuously decreasing temperature with conver¬ 
sion cause the strong decline of reaction rate towards the gasifier 
outlet. The gas outlet concentrations are: 22.1% H 2 , 50.8% CO, 
13.4% H 2 0, 9.1% C0 2 and 4.6% N 2 . These are in the expected range 
for entrained flow gasification of lignite [41 ]. 

The GSP gasifier for a thermal input of 500 MW has a length of 
5.25 m and a diameter of 2.9 m [41 ]. The gasifier is designed for 
various fuels including bituminous coals that are less reactive. 
Therefore a downsizing of the gasifier is expected when fuels with 
a higher reactivity are used. The simulation proposes a smaller gas¬ 
ifier size for lignite that has a high reactivity. This is an indicator 
that the gasification model gives meaningful results when trans¬ 
ferred to an industrial scale. 


Parameter 


Thermal input 
Pressure 
O/C ratio 
Oxygen purity 
Inlet temperature 
Carrier gas (N 2 ) 


500 MW 

2.5 MPa [42] 

0.95 

99% [43] 

90 °C 

0.07 m 3 /kg fuel [43] 


4.3. Optimisation of the gasifier operation 

The gasifier model is used to analyse the influence of operation 
parameters and boundary conditions on the gasifier performance. 
The gasifier operation aims at an optimisation of cold gas efficiency 
and fuel conversion. The effects of stoichiometry, of steam or C0 2 
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addition to the burner, of a variation in carrier gas flow rate, of fuel 
particle size, and of heat losses are evaluated in a sensitivity anal¬ 
ysis. The oxygen to fuel stoichiometry is the most crucial parame¬ 
ter. Its influence on CGE and fuel conversion is shown in Fig. 8. 

There is an optimum O/C ratio around 0.88, and the CGE drops 
rapidly at both sides of this value. However, an increase in oxygen 
addition improves the fuel conversion due to the higher tempera¬ 
ture in the gasifier. If the complete conversion is required or heat 
losses occur, the practical O/C ratio is higher than the CGE opti¬ 
mum. Furthermore, the gasifier outlet temperature is decreased 
by the reduction of the stoichiometry. At an O/C ratio of 0.88 the 
outlet temperature is only 1215 °C which is probably too low to 
achieve a steady slag flow. If more information about the slag vis¬ 
cosity at a certain temperature is available, the O/C ratio can also 
be adapted to the slag requirements. 

A further important process parameter is the fuel particle size. 
In entrained flow gasification the particle size has to be small en¬ 
ough to reduce mass transport limitations and achieve fuel conver¬ 
sion in the given time interval. However, a small particle size 
results in a higher energy consumption of the fuel mill. It is helpful 
to analyse the influence of particle size on the gasifier perfor¬ 
mance. This is shown in Fig. 9. 

Both CGE and fuel conversion drop with increasing particle size. 
Up to a particle size of about 200 pm the effect is small, but larger 
particles lead to a massive increase in mass transport limitations 
within the pore structure and the boundary layer. The average par¬ 
ticle effectiveness factor is an indicator of the severity of the mass 
transport limitations within the char particle structure. It is 0.82 
for 60 pm particles, 0.63 for 100 pm particles, but decreases down 
to 0.02 for 500 pm particles. Therefore the surface area towards the 
centre of larger particles does not contribute to char conversion, 
and fuel conversion is significantly decreased. 

The influence of particle size is in agreement with particle 
diameter requirements given by Schingnitz and Mehlhose [42], 


For hard coal with high reactivity and hard brown coal the GSP 
process requires that minimum 94% of all particles are smaller than 
250 pm to achieve a conversion of 99%. This requirement is also 
suggested by the simulation data in Fig. 9. 

4.4. Analysis of the reaction regime in entrained flow gasification 

Experimental data that are found in the literature are limited by 
an upper temperature of 1500 °C. Only a handful of studies is found 
that discus the reaction regime of char gasification under entrained 
flow conditions. Hodge [16] observe the transition from Regime I to 
Regime II at or around 1000-1100 °C and assume that in entrained 
flow gasifiers at a temperature above 1200 °C the char-C0 2 reac¬ 
tion will occur under Regime II conditions. In other experimental 
studies [11-13] the transition temperatures to Regime II condi¬ 
tions are between 1000 °C and 1400 °C for the C0 2 and H 2 0 gasifi¬ 
cation and a transition to Regime III is not considered. In solid fuel 
combustion, the transition to Regime III conditions occurs at lower 
temperature than in entrained flow gasification due to the very fast 
reaction rate of the char-0 2 reaction. In entrained flow gasifiers a 
transition to Regime III is expected at very high temperature. 

The experimental data up to 1600 °C that are shown here can¬ 
not be directly used to determine the reaction regime. However, 
the model framework derived on the basis of these data can be ap¬ 
plied to analyse the reaction regime. 

The gas concentrations in the bulk gas phase and at the char 
particle outer surface are an indicator of the reaction regime. These 
gas concentrations for the char-H 2 0 reaction within the 500 MW 
gasifier for coal R are shown in Fig. 10. Within the first metre of 
the gasifier length there is a large difference between the H 2 0 bulk 
and surface concentrations. This is due to mass transfer limitations 
within the boundary layer of the char particles and is a clear indi¬ 
cator of Regime III conditions. After a reaction length of 2 m there 
is no H 2 0 concentration gradient in the boundary layer. The mass 



Fig. 8. Influence of oxygen to fuel stoichiometry (O/C ratio, 0% value is 0.95) on: (a) cold gas efficiency and (b) fuel conversion of a 500 MW entrained flow gasifier. 
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transport limitation within the pore structure is expressed by the 
effectiveness factor that is also shown in Fig. 10. At a high temper¬ 
ature in the initial stages of reaction the effectiveness factor is 0.02 
indicating a strong limitation by pore diffusion. Therefore the reac¬ 
tion rate is initially limited by boundary layer and pore diffusion. 
Then the effectiveness factor continuously increases and ap¬ 
proaches almost unity at the gasifier outlet. Therefore, the reaction 
rate approaches Regime I conditions in the later stages of 
conversion. 

The decrease in reaction rate and the continuous transfer from 
Regime III to Regime 1 conditions is caused by the combination of 
thermal deactivation, loss of surface area and the continuously 
decreasing temperature with conversion. 

The temperature profile and the local reaction rate in the gas¬ 
ifier can be used to derive a diagram similar to the Arrhenius plot. 
The logarithm of the apparent reaction rate [g/(gs)] is plotted ver¬ 
sus the inverse temperature in Fig. 11. In the classical Arrhenius 
plot, only the influence of temperature is shown whereas all other 
parameters are kept constant. Fig. 11 is directly derived from the 
simulation and the reaction rate in the diagram is also influenced 
by variations of other parameters (gas concentrations, surface area, 
char deactivation, etc.). But due to its exponential influence the 
temperature is the most important parameter. Due to the influence 
of other parameters the activation energy of the reaction cannot be 
derived from the diagram. The slope of the curve (observed 
activation energy) in Fig. 11 can be used to discuss the impact of 
mass transfer limitation on char conversion. 



low temperature is indicated by the observed activation energy). 


At low temperatures near the gasifier outlet the temperature 
influence is described by an observed activation energy of 239 kj/ 
mol. This is slightly above the activation energy for the intrinsic 
char-C0 2 and char-H 2 0 reactions because further effects (surface 
area, char deactivation, etc.) also influence the observed reaction 
rate. The Regime I conditions near the gasifier outlet at low tem¬ 
perature are confirmed due to the steep decline of the curve (high 
observed activation energy). 

The slope of the curve continuously decreases at higher 
temperatures, but a defined transition temperature to Regime II 
conditions is not observed. In ideal Regime II conditions the ob¬ 
served activation energy is half the activation energy under Regime 
I conditions. This is the case in the medium region of the gasifier. At 
the gasifier top, the maximum temperature during char conversion 
occurs and the observed activation energy is only 62 kj/mol. As this 
value is very low compared to the observed activation energy un¬ 
der Regime I conditions, the reaction seems to occur in Regime III. 
The plot of apparent rate versus inverse temperature confirms the 
mass transport limitations indicated by the variations of gas con¬ 
centrations shown in Fig. 10. 

The particle size is another important parameter that deter¬ 
mines the reaction regime. At a larger particle size the reaction 
shifts further to Regime III conditions. In the evaluation so far a 
particle size of 100 pm has been used. If the particle diameter is in¬ 
creased to 200 pm the initial activation energy is 49 kj/mol, and at 
a particle size of 500 pm it is only 26 kj/mol, indicating Regime III 
conditions. 


5. Conclusion 

A gasification model is presented that allows the simulation of 
fuel reaction behaviour under entrained flow conditions, the de¬ 
sign of industrial scale gasifiers, and the investigation of the influ¬ 
ence of their operating conditions. 

Lignite is used in the investigations and the model is validated 
using entrained flow data at high pressure and high temperature. 
A high conversion of the coal is achieved in the experiment. The 
transfer of the data to the industrial scale suggests that adequate 
conversion can be achieved under entrained flow conditions using 
a compact gasifier design for a thermal input of 500 MW. Therefore 
the experiments and the simulations indicate that the lignite is a 
suitable fuel for entrained flow gasification. 

In the simulation the operation condition of the gasifier are var¬ 
ied and an optimum stoichiometry is found. If less oxygen is fed to 
the burner the fuel does not achieve complete conversion. At a 
higher stoichiometry the percentage of fuel that is combusted is in¬ 
creased and the cold gas efficiency of the process is reduced. 

The simulation shows that all three reaction regimes are impor¬ 
tant during entrained flow gasification. In order to predict reaction 
behaviour under entrained flow conditions, evaluation of the ef¬ 
fects of operation parameters and boundary conditions in all three 
regimes is required and intrinsic reactivities as well as mass trans¬ 
port limitations must be included. 

A simple one-dimensional plug flow reactor is chosen to simu¬ 
late a larger scale entrained flow gasifier. The temperature and gas 
profiles within an entrained flow gasifier are shown. However, de¬ 
tails like the burner zone cannot be addressed in detail. The model 
structure is appropriate for the implementation in a CFD software 
package. Further work will apply the gasification model in 
combination with numerical methods to show and improve the 
operation of entrained flow gasifiers. 

Further investigations should also focus on the derivation of the 
model parameters and the validation of the model for other fuels. 
Due to its general approach the gasification model is valid for all 
solid fuels ranging from high rank coal to biomass. 
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